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We describe the primeval inflationary phase of the early Universe within a quantum field theo¬ 
retical (QFT) framework that can be viewed as the effective action of vacuum decay in the early 
times. Interestingly enough, the model accounts for the “graceful exit” of the inflationary phase 
into the standard radiation regime. The underlying QFT framework considered here is Super¬ 
gravity (SUGRA), more specifically an existing formulation in which the Starobinsky-type inflation 
(de-Sitter background) emerges from the quantum corrections to the effective action after integrating 
out the gravitino fields in their (dynamically induced) massive phase. We also demonstrate that the 
structure of the effective action in this model is consistent with the generic idea of renormalization 
group (RG) running of the cosmological parameters, specifically it follows from the corresponding 
RG equation for the vacuum energy density as a function of the Hubble rate, pa{H). Overall our 
combined approach amounts to a concrete-model realization of inflation triggered by vacuum decay 
in a fundamental physics context which, as it turns out, can also be extended for the remaining 
epochs of the cosmological evolution until the current dark energy era. 
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I. INTRODUCTION 

In the last two years we have witnessed extraordi¬ 
nary developments on experimental tests of inflation¬ 
ary models [l|, based on studies of photons in the cos¬ 
mic microwave background radiation. In particular, the 
results of Planck collaboration 0 and the associated 
non-observation of B-mode polarizations of primordial 
light fluctuations, have imposed very stringent restric¬ 
tions on single scalar-field models of slow-roll inflation, 
allowing basically models with very low tensor-to-scalar 
fluctuation ratio r = n'r/ug <C I, with a scalar spec¬ 
tral index rig ~ 0.96 and no appreciable runnin g. In 
fact the upper bound set by Planck Collaboration on 
this ratio, as a consequence of the non-observation of 
B-modes, is r < O.II, but their favored regions point 
towards r < 10“^. This is a feature that character¬ 
izes the so-called Starobinsky-type (or i?^-inflation, with 
R denoting the scalar space-time curvature) inflationary 
models [^. The estimated energy scale Ej of inflation, 
which in inflaton-type models is related to the - ap¬ 
proximately constant- scalar potential during inflation 
through Ej = reads [l|: 

Ej = c 2.1 X X (—) GeV , 

( 1 ) 

where Mpi = I/VSttG ~ 2.43 x 10^® GeV is the reduced 
Planck mass (G being the Newtonian constant). The 
upper bound r < 0.11 placed by the Planck Collaboration 


implies 

Hi = 1.05 X 10i'‘ GeV < 0.78 x 10i'‘ GeV . 

( 2 ) 

The above can be rephrased as Hj/mp ^ 6.39 x 10 
where mp = IjG^I'^ = ~ 1.22 x 10^® GeV is the 

Planck mass in natural units. This result is consistent 
with the well known CMB bound Hj/mp < 10“^ on the 
temperature fluctuations induced by the tensor modes. 
As we will see, the actual value of H during inflation for 
the class of models under study satisfies H < Hj, and 
hence the GMB bound is preserved by them. 

The recent joint BICEP2-Planck analysis @ confirmed 
the early Planck result, namely the likelihood curve for 
r yields an upper limit r < 0.12 at 95%. Moreover, the 
present BIGEP2-Planck data are consistent with a scalar 
spectral index rig ~ 0.96 and no appreciable running, in 
agreement with the previous Planck data 0. Using the 
aforementioned new upper limit rmax = 0.12, the Hubble 
parameter during slow-roll inflation Hj is estimated to 
be below 

^Bicep2+Planck ^ q ^ ^q14 ^ 

and hence Hj/mp ^ 6.64 x 10“®. Because of the low 
significance of the new limit on r, the possibility that 
r is actually much smaller than the current upper limit 
^■max remains as natural as it was before. In fact, nothing 
actually prevents at present that the typical value of the 
tensor to scalar ratio can be, for example, r = 0(10“®), 
and in this sense the Starobinsky-type scenarios can still 
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be considered as a serious possibility to describe the in¬ 
flationary universe. Following this point of view, we con¬ 
tinue in this paper with the investigation of Starobinsky- 
like models as potential candidates for realistic imple¬ 
mentation of inflation compatible with the data. 

In previous publications one of us (N.E.M.) with col¬ 
laborators 0 discussed the dynamical breaking of su¬ 
pergravity (SUGRA) theories via gravitino condensation, 
and demonstrated Q the compatibility of this scenario 
with Starobinsky-like inflationary scenarios. As we 
discussed, this phase is characterized by the dynami¬ 
cal emergence of a de-Sitter background. As argued in 
the Starobinsky-type inflation appears much more 
natural (from the point of view of the order of the pa¬ 
rameters involved) than a hill-top inflation scenario @ 
in which the gravitino condensate itself is the inflaton 
field. In the latter, very large values of the wave function 
renormalization of the condensate field are required to 
ensure slow-roll inflation if one insists on (phenomeno¬ 
logically realistic) sub-Planckian supersymmetry break¬ 
ing scales. It is important to notice at this point that 
in the original Starobinsky model Q the E? terms cru¬ 
cial for inflation arise from the conformal anomaly in the 
path integral of massless (conformal) matter in a de Sit¬ 
ter background, and thus their coefficient is arbitrary, 
and can only be fixed phenomenologically. A similar, al¬ 
though not identical, situation occurs in the context of 
anomaly-induced inflation , where the term Ef is 

absent at the classical level but is generated from the con¬ 
formal anomaly. In this case, however, the coefficient of 
□i? (entering the /3-functions and controlling the stability 
of inflation) presents also some arbitrariness which can 
only be fixed by a special renormalization condition. Par 
contrast, in the considered SUGRA scenario, such terms 
arise in the one-loop effective action of the gravitino con¬ 
densate field, evaluated in a de Sitter background, after 
integrating out massive gravitino fields, whose mass was 
generated dynamically. The order of the de Sitter cosmo¬ 
logical constant, A > 0 that breaks supersymmetry, and 
the gravitino mass are all evaluated dynamically (self- 
consistently) in our approach from the minimization of 
the effective potential. Thus, the resulting E^ coefficient, 
which determines the phenomenology of the inflationary 
phase, is calculable Q- 

Also very important for our considerations is the 
framework of the running vacuum model (RVM) [lol - [l3l| 
- see MM and references therein for a recent detailed 
exposition. The implications of these dynamical vacuum 
models have recently been analyzed both for the early 
universe as well as for the phenomenology of the 

current universe [2ll - l^ - see also [23 - 1^ for previous 
analyses. 

In regard to the early universe we emphasize that the 
RVM defines a class of non-singular inflationary scenar¬ 
ios with graceful exit into the standard radiation regime. 
These models are related to Starobinsky inflation mod¬ 
els, although they are not equivalent. We will discuss 


in this paper the correspondence between them, and 
most particularly with the dynamically broken SUGRA 
model with gravitino condensation that we have men¬ 
tioned above. It is especially remarkable that such spe¬ 
cific implementation of the SUGRA model leads, as we 
will show in this paper, to the effective behavior of the 
RVM with calculable coefficients. In this way the former 
automatically benefits from the successful consequences 
of the latter. Let us mention that the RVM also pro¬ 
vides some important clues for alleviating the cosmolog¬ 
ical constant problem M- 

Finally, we should like to mention that the RVM’s 
have been tested against the wealth of accurate 
SNIa-bBAO-fiL(z)-fLSS-bBBN-fGMB data - see [H for 
a recent summary review - and they turn out to provide 
a quality fit that is significantly better than the AGDM. 
This fact has become especially prominent in the light of 
the most recent works [sll, [s^. Therefore, there is every 
motivation for further investigating these dynamical vac¬ 
uum models from different perspectives, with the hope 
of finding possible connections with fundamental aspects 
of the cosmic evolution. In point of fact, this is the main 
aim of this work. 

The structure of the article is as follows. The gen¬ 
eral framework of the RVM is introduced in In] The basic 
theoretical elements of the Starobinsky inflation are pre¬ 
sented in section cni The main properties of the dynam¬ 
ical breaking of local SUGRA theory and its connection 
to Starobinsky-type inflation are reviewed in section HVl 
In section |V] we demonstrate how the RVM describes the 
effective framework of the Starobinsky Q and the dy¬ 
namically broken SUGRA Q models at the inflationary 
epoch. Finally, our conclusions are summarized in sec¬ 
tion |VT1 


II. RUNNING VACUUM: A NATURAL ARENA 
FOR VACUUM DECAY IN COSMOLOGY 

It is the purpose of this work to go one step fur¬ 
ther from demonstrating compatibility of the dynami¬ 
cally broken SUGRA scenario @,0 with inflation and 
discuss the possibility of a dynamical evolutionof the 
inflationary phase ground state to the standard radia¬ 
tion regime within the context of the running vacuum 
model (RVM) of the cosmic evolution utilizing an effec¬ 
tive “renormalization group (RG) approach”, see [13, [HI 
- and 0111 for comprehensive expositions. Specifi¬ 
cally, we wish to show that the behavior of the afore¬ 
mentioned SUGRA scenario effectively mimics the RVM. 
Once this link is elucidated, the general “decaying” 
vacuum description inherent to the RVM formulation 
allows to smoothly connect inflation to the standard 
Fridman-Lemaitre-Robertson-Walker (FLRW) radiation 
era, which subsequently proceeds into a matter and dark 
energy domination in the present era, in which it still 
carries a mild dynamical behavior compatible with the 
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current cosmological data MM- Such an expansion 
history of the Universe has been put forward in previ¬ 
ous works by the authors in various collaborations and 
contexts, e.g. non-equilibrium string-inspired cosmolo¬ 
gies [13 or conventional field-theoretic cosmologies in 
which the above mentioned RVM is extensively applied 
for the study of the early cosmic history [IGI - IiM . 

In the effective RG approach underlying the RVM one 
can write down an evolution equation for the effective 
vacuum energy density pA(i) = PAifJ-cit)), treated as a 
dynamical quantity whose cosmic time evolution is in¬ 
herited from its dependence on a characteristic cosmic 
scale variable Pc = Pc(i)- This variable plays the role 
of running (mass) scale of the renormalization group ap¬ 
proach, and a natural candidate for such scale in FLRW 
cosmology is the Hubble parameter H{t). Therefore the 
proposed RG equation is [ij]: 


dp\{t) _ 
dlni/2 “ ■(47r) 


.E 




W6 

iMf H'^ + biH'^ + Ci 


( 4 ) 


In general /i^ can be associated to a linear combination 
of and H and the variety of terms appearing on the 
r.h.s. of (m can be richer [l^ ■ but the canonical possi¬ 
bility is the previous one and hereafter we restrict to it. 
The coefficients Oj, ... appearing in (|3]) are dimen¬ 

sionless and receive contributions from loop corrections 
of boson and fermion matter fields with different masses 
Mi. It must be stressed that the general covariance of 
the action US [mil necessitates the appearance of only 
even powers of the (cosmic-time t dependent) Hubble pa¬ 
rameter H{t) on the right-hand-side of (HJ. For a specific 
framework where the above RG is concretely realized and 
the /3-function coefficients can be computed, see M- 

We note at this stage that, if the evolution of the Uni¬ 
verse is restricted to eras below the Grand Unified The¬ 
ory (GUT) scale, then for all practical purposes it is at 
most the terms (those with dimensionless coefficients 
hi) that can contribute significantly. The term is of 
course negligible at this point, and the higher powers of 
iJ" for n = 6,8,.. are suppressed by the corresponding 
inverse powers of the heavy masses Mi, which go to the 
denominator, as required by the decoupling theorem. In 
the scenarios of dynamical breaking of local supergrav¬ 
ity discussed in ref. @-H, the breaking and the associ¬ 
ated inflationary scenarios could occur around the GUT 
scale, in agreement with the inflationary phenomenol¬ 
ogy suggested by the Planck satellite data d, provided 
Jordan-frame supergravity models (with broken confor¬ 
mal symmetry) are used, in which the conformal frame 
function acquired, via appropriate dynamics, some non 
trivial vacuum expectation value. For these situations, 
therefore, corrections in (j4]) involving higher powers than 
will be ignored. 

In the next sections, after revising the general frame¬ 
work of Starobinsky inflation, we shall compute p\ in 
such supergravity models and study their evolution from 
the exit from the Starobinsky inflationary phase, that 


occurs in the massive gravitino phase, until today. The 
computation of pA will be made via the corresponding 
calculation of the one-loop effective action after massive 
gravitinos are integrated out in a path integral. Then, 
an identification of the effective equation of state can be 
derived by integrating ®, following the approach of 0- 
IH, [11. Before doing so, it is instructive to review first 
the emergence of Starobinsky-type inflation. 


III. GENERIC STAROBINSKY INFLATION 

Starobinsky inflation is the oldest model of inflation Q , 
prior to the traditional, scalar-field-based, inflaton mod¬ 
els. It is characterized for being able to realize the de 
Sitter (inflationary) phase from the gravitational field 
equations derived from a four-dimensional action that 
includes higher curvature terms, specifically of the Wpe 
involving the quadratic curvature correction ^ Hu: 

J (i? + /3i?") , ^ = 

Here = SttG = 1/Mpj (in units of /i = c = 1 we 
are working on), G = l/rup is Newton’s (gravitational) 
constant in four space-time dimensions, with mp the 
Planck mass, and AI is a constant of mass dimension 
one, characteristic of the model. Notice that the curva¬ 
ture terms in the action are just the dimension-4 combi¬ 
nation m|.i?/167r-|-i?^mp/(6AI^). With this normaliza¬ 
tion, A4 gives the value of the so-called scalaron mass. 
The smaller is Ai in Planck mass units (i.e. the larger 
is the dimensionless parameter rup/AI^ in front of R^) 
the longer is the inflationary time (cf. Fig. 2 of ill)- 
Of course Ai cannot be much below the natural scale 
of inflation, and in fact it should be of the same order, 
i.e. A4 ^ Mx, where Mx is some GUT scale below the 
Planck mass. Typically Mx ^ 10^® GeV ^ 10“^mp. 

The most relevant feature of this model is that infla¬ 
tionary dynamics is driven by the purely gravitational 
sector, through the R^ terms. From a microscopic point 
of view, these terms can be viewed as the result of quan¬ 
tum fluctuations (at one-loop level) of conformal (mass¬ 
less or high energy) matter fields of various spins, which 
have been integrated out in the relevant path integral in a 
curved background space-time [11 . The model in fact is 
to be understood in the context of QFT in curved space- 
time. The quantum mechanics of this model, by means 
of tunneling of the Universe from a state of “nothing” 
to the inflationary phase of ref. [1 has been discussed 
in detail in [11. The above considerations necessitate 
truncation to one-loop quantum order and to curvature- 
square (four-derivative) terms, which implies that there 


^ Our metric signature is -t,-t) and the definitions of the 

Ricci and Riemann curvature tensors are and 

R^vp — dv — .. ■, respectively, i.e. we follow the exact three- 
sign conventions (-l-,-|-,+) of Misner-Thorn-Wheeler |34| . 
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must be a region of validity for curvature invariants such 
that 0(i?^/TOp) <C 1. Recalling that R ^ 12 Hf in the 
inflationary phase (where Hi is the nearly constant Hub¬ 
ble rate in that phase), we observe that this is indeed 
a condition satisfied in phenomenologically realistic sce¬ 
narios of inflation 0,i, for which the inflationary Hub¬ 
ble scale Hi is typically constrained to obey (H)) (Planck 
data 0) or ([3|) (BICEP2 data 0), which are at present 
essentially the same. 

Although the inflation in this model is not driven 
by fundamental rolling scalar fields, nevertheless the 
model ([S]) (and for that matter, any other model where 
the Einstein-Hilbert space-time Lagrangian density is re¬ 
placed by an arbitrary function f{R) of the scalar cur¬ 
vature) is conformally equivalent to that of an ordinary 
Einstein-gravity coupled to a scalar held with a potential 
that drives inflation [^. To see this, one firstly linearises 
the R^ terms in ([5]) by means of an auxiliary (Lagrange- 
multiplier) held a{x), before rescaling the metric by a 
conformal transformation and redefining the scalar field 
(so that the final theory acquires canonically-normalised 
Einstein and scalar-field terms): 

gu’' g^’' = (1 + 2^d(a:)) , 

a (a;) -)■ Kip{x) = hi{l + 2j5a (a;)) , (6) 

where again k = -v/SttG. These steps may be understood 
schematically via 

J d*x^ (R + PR^) ( 7 ) 

^ J d‘^X\/—g ^(1 -f 2 /3 d (a:)) R — P a[x)^'^ 

^ j d^x^/^ {R^ , 


where the arrows have the meaning that the correspond¬ 
ing actions appear in the appropriate path integrals. The 
ensuing effective potential Ves(j-p) is given by: 


VM = 


3M^ 1-e 


— e vT'' 


4k^ 


( 8 ) 


One can check that the mass of the scalaron, which can 
be seenn as the new gravitational degree of freedom that 
the conformal transformation was able to elucidate from 
the Starobinsky action, is indeed given by the parameter 
A4 = ^Sn/SP: 


d'^Vesjip) 

dip^ 


= M^. 

V ?—0 


(9) 


Note that for (p = 0 one has 5 = 0 and the two con¬ 
formally equivalent metrics coincide at this point. The 
effective potential for the scalar d.o.f. that conformally 
replaces the effect of the R? term is plotted in Fig. [T] We 
observe that V{ip) is sufficiently flat for Kp> = (p/Mpi ^ 1 





FIG. 1. The effective potential ([8]) of the collective scalar field 
p that describes the one-loop quantum fluctuations of matter 
fields, leading to the higher-order scalar curvature corrections 
in the Starobinsky model for inflation ([5]). Notice that accord¬ 
ing to (|8| we have P = The potential is sufBciently 

flat for np 1 to ensure slow-roll conditions for inflation are 
satisfied, in agreement with the Planck data, for appropriate 
values of the scale 1 /oc M (which sets the overall scale of 
inflation in the model). 


(i.e. for sufficiently large values of p as compared to the 
reduced Planck scale) to produce phenomenologically ac¬ 
ceptable inflation. Obviously the scalaron field p is ef¬ 
fectively playing the role of the inflaton in this context. 
The difference with the usual inflaton is that p is not 
a new scalar d.o.f. imported from outside the gravita¬ 
tional action, but just an integral part of it, namely, it 
is just a gravitational d.o.f. that describes in an effec¬ 
tive (and very convenient) way the ^ Rp term of (|5|) . 
The Starobinsky model based on the action dS]) indeed 
fits excellently the Planck data on inflation Q, and also 
the corresponding data from the joint BICEP2-Planck 
analysis Q . 

Quantum-gravity corrections in the original Starobin¬ 
sky model ([5]) have been considered recently in from 
the point of view of an exact renormalisation-group anal¬ 
ysis [1^. It was shown that the non-perturbative beta- 
functions for the ‘running’ of Newton’s ‘constant’ G and 
the dimensionless inverse R^ coupling Kpp~^ ^ A4^/M|>j 
in m imply an asymptotically safe Ultraviolet (UV) 
fixed point for the former (that is, G{k —>■ oo) —>■ con¬ 
stant, for some 4-momentum cutoff scale fc), in the spirit 
of Weinberg [i^, and an attractive asymptotically-free 
{K^P~^{k —>■ oo) —>■ 0) point for the latter. In this sense, 
the smallness of the (inverse) R^ coupling, required for 
agreement with inflationary observables [^, is naturally 
ensured by the presence of the asymptotically free UV 
fixed point. 

The agreement of the model of Q with the Planck data 
triggered an enormous interest in the current literature, 
and indeed Starobinsky inflation has been revisited from 
various points of view, such as its connection with no- 
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scale supergravity [4l| and (super)conformal versions of 
supergravity and related areas In the latter works 
however the Starobinsky scalaron field is fundamental, 
arising from the appropriate scalar component of some 
chiral superfield that appears in the superpotentials of 
the model. 

Although of great value, illuminating a strong con¬ 
nection between supergravity models and inflationary 
physics, and especially for explaining the low-scale of in¬ 
flation compared to the Planck scale, these works con¬ 
tradict the original spirit of the Starobinsky model ([S]) 
where, as mentioned previously, the higher curvature 
corrections are viewed as arising from quantum fluctua¬ 
tions of matter fields in a curved space-time background 
such that inflation is driven by the pure gravity sector in 
the absence of fundamental scalars. On the other hand, 
the scenario of Q, in which a Starobinsky-type inflation 
arises in the massive gravitino phase of SUGRA models, 
after integrating out the massive degrees of freedom is 
in the same spirit of Starobinsky, and even better, in the 
sense that the model does not have to assume dominance 
of conformal matter during inflation. 

We next proceed to summarize the construction of 
the one-loop effective action of the massless degrees 
of freedom after massive gravitino integration in this 
dynamically-broken SUGRA model with spontaneous 
breaking of global supersymmetry (SUSY) 

IV. STAROBINSKY-TYPE INFLATION IN 
DYNAMICALLY BROKEN SUGRA 

Dynamical breaking of SUGRA, in the sense of the 
generation of a mass for the gravitino field , whilst the 
gravitons remain massless, occurs in the model as a re¬ 
sult of the four-gravitino interactions characterizing the 
SUGRA action, arising from the torsionful contributions 
of the spin connection, characteristic of local supersym¬ 
metric theories. 

Our starting point is the J\f = 1 D = A (on-shell) action 
for ‘minimal’ Poincare supergravity in the second order 
formalism (d^ : 

SSG = J d'^Xe (e) - + Aorsion^ , 

( 10 ) 

k^=87tG, = , 7-P = i[7-,7P] , 

where R{e) and D^tpp = duijjp -I- jUJ^ab (e) 7“^V’p are de¬ 
fined via the torsion-free connection; and, given the gauge 
condition 7 -^ = 0 , 

2 

Aorsion = -y > 

( 11 ) 

arising from the fermionic torsion parts of the spin con¬ 
nection. Extending the action off-shell requires the addi¬ 


tion of auxiliary fields to balance the graviton and grav¬ 
itino degrees of freedom. These fields however are non¬ 
propagating and may only contribute through the devel¬ 
opment of scalar vacuum expectation values, which would 
ultimately be resummed into the cosmological constant. 

Making further use of the above gauge condition to¬ 
gether with the Fierz identities (as detailed in @), we 
may write 

^” 2 *°" = + Aps a) + Apv (^U^7®7 pUp) 

( 12 ) 

where the couplings As, Aps and Apv express the free¬ 
dom we have to rewrite each quadrilinear in terms of the 
others via Fierz transformation. This freedom in turn 
leads to a known ambiguity in the context of (perturba¬ 
tive) mean field theory and can only be resolved by 
a non-perturbative treatment. 

Specifically, we wish to linearise these four-fermion in¬ 
teractions via suitable auxiliary fields, e.g. 

^k'^Xs (^'tp^'lpp^ - cr K (^^V’p) - , (13) 

where the equivalence (at the level of the action) fol¬ 
lows as a consequence of the subsequent Euler-Lagrange 
equation for the auxiliary scalar a. Our task is then to 
look for a non-zero vacuum expectation value (tr) which 
would induce as an effective mass 7713/2 ^ a k a/As for the 
gravitino. This is however complicated by the fact that 
our coupling As into this particular channel is, by virtue 
of Fierz transformations, ambiguous at a perturbative 
level and, as mentioned, in order to fix them a fully non- 
perturbative treatment of SUGRA-like models would be 
required, which are not currently at hand. Nevertheless, 
there is another way out [1, whereby the Fierz am¬ 
biguities may be absorbed by dilaton-expectation-value 
shifts in an extension of A/" = 1 SUGRA which incorpo¬ 
rates local supersymmetry in the Jordan frame, enabled 
by an associated dilaton superfield [d^. The (logarithm 
of the) scalar component ip of the latter can be either a 
fundamental space-time scalar mode of the gravitational 
multiple!, i.e. the trace of the graviton (as happens, for 
instance, in supergravity models that appear in the low- 
energy limit of string theories), or a composite scalar field 
constructed out of matter multiplets. In the latter case 
these could include the standard model fields and their 
superpartners that characterise the Next-to-Minimal Su¬ 
persymmetric Standard Model [d^ , which can be consis¬ 
tently incorporated in such Jordan frame extensions of 
SUGRA. 

Upon appropriate breaking of conformal symmetry, in¬ 
duced by specific dilaton potentials (which we do not 
discuss here), one may then assume that the dilaton field 
acquires a non-trivial vacuum expectation value ((/>) yf 0, 
thus absorbing any ambiguities in the value of the ap¬ 
propriate coefficient Ag induced by Fierz ca). One con¬ 
sequence of this is then that in the broken conformal 
symmetry phase, the resulting supergravity sector, upon 
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passing (via appropriate field redefinitions) to the Ein¬ 
stein frame is described by an action of the form 
but with the coupling of the gravitino four-fermion inter¬ 
action terms being replaced by 

= As , (14) 

while the Einstein term in the action carries the stan¬ 
dard gravitational coupling 1/2^^. Eor phenomenologi¬ 
cal reasons, associated with gravitino masses in the ball¬ 
park of Grand Unified Theory (GUT) scales, one must 
have k ^ K. This is assumed to be guaranteed by ap¬ 
propriate microscopic dilaton potentials that break the 
(super) conformal symmetry of the Jordan-frame SUGRA 
appropriately. 

To induce the super-Higgs effect [l^ we couple to the 
action m the Goldstino associated to global supersym¬ 
metry breaking via the addition of 

Cx = f det (= f + ^ 

(15) 

where A is the Goldstino, a// expresses the scale of global 
supersymmetry breaking, and ... represents higher or¬ 
der terms which may be neglected in our weak-field ex¬ 
pansion of the determinant. It is worth emphasising at 
this point the universality of (HID; any model containing 
a Goldstino may be related to [L\ via a non-linear trans¬ 
formation [ 4 ^, and thus the generality of our approach 
is preserved. 

Upon the aforementioned gauge choice for the grav¬ 
itino field = 0 > appropriate redefinition, 

one may eliminate any presence of the Goldstino field 
from the final effective action describing the dynamical 
breaking of local supersymmetry, except the cosmologi¬ 
cal constant term in (11511 . which serves as a reminder 
of the pertinent scale of supersymmetry breaking. The 
non-trivial energy scale this introduces, along with the 
disappearance (through field redefinitions) of the Gold¬ 
stino field from the physical spectrum and the concomi¬ 
tant development of a gravitino mass, characterises the 
super-Higgs effect. 

The linearisation of the four-gravitino terms (HI, 
when combined with the term of the super-Higgs effect 
implies a tree-level cosmological constant 

( 16 ) 

which must be negative due to the incompatibility of su¬ 
pergravity with de Sitter vacua (notice that in our con¬ 
ventions both a and / have dimension +2 in natural 
units). 

The one-loop effective potential for the scalar gravitino 
condensate field a{x) (with vacuum expectation value 
tJc oc tias a double-well shape as a function of 

a{x) which is symmetric about the origin (c/. fig. [5]), as 
dictated by the fact that the sign of a fermion mass does 


not have physical significance. Dynamical generation of 
the gravitino mass occurs at the non-trivial minima cor¬ 
responding to (Tc 0. The potential of the CTc field is 
also flat near the origin, and this has been identified in 
Q with a first inflationary phase. 

In the one-loop effective potential was derived by 
first formulating the theory on a curved de Sitter back¬ 
ground [ 53 , HH , with cosmological constant (one-loop in¬ 
duced) A > 0, not to be confused with the (negative) 
tree-level one Aq dni), and then integrating out spin-2 
(graviton) and spin 3/2 (gravitino) quantum fluctuations 
in a given class of gauges (physical), before considering 
the flat limit A —>■ 0 in a self-consistent way. The de¬ 
tailed analysis in performed in the physical gauge, 
has demonstrated that the dynamically broken phase is 
then stable (in the sense of the effective action not being 
characterized by imaginary parts) provided the scale of 
the gravitino condensate is equal or below the scale of 
spontaneous breaking of global SUSY: 

(17) 

which guarantees the aforementioned result on the neces¬ 
sity of the negative nature of the tree-level cosmological 
constant mi). 

The former result demonstrates the importance of the 
existence of global SUSY breaking scale for the stabil¬ 
ity of the phase where dynamical generation of gravitino 
masses occurs, which was not considered in the previous 
literature [i^. In super-conformal versions of SUGRA, 
e.g. those in ref. [^, 1^, phenomenologically realistic 
scales for and gravitino mass of order of the GUT 
scale, appear for appropriate values of the expectation 
value of the conformal factor. These imply inflationary 
scenarios in perfect agreement with the Planck data @,3, 
on equal footing to the original Starobinsky model. 

In 0 we considered an extension of the analysis of @ 
to the case where the de Sitter parameter A is pertur- 
batively small compared to Wp, but non zero, so that 
truncation of the series to order A^ suffices. This is in 
the spirit of the original Starobinsky model Q , with the 
role of matter fulfilled by the now-massive gravitino field. 
Specifically, we were interested in the behavior of the 
effective potential near the non-trivial minimum, where 
cr ~ (Tc is a non-zero constant (c/ Eig. [2]). The one-loop 
effective potential, obtained by integrating out [13, Ull 
gravitons and (massive) gravitino fields in the scalar 
channel (after appropriate euclideanisation), may be ex¬ 
pressed as a power series in A: 

r ~ Sci + (“f + ^ 

+ {pi2 + 012 ) A^ -I-...) , (18) 

where ^cl denotes the classical action with tree-level cos¬ 
mological constant Aq (to be contrasted with the one- 
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FIG. 2. Generic shape of the one-loop effective potential 
(in dimensionless units Vefi , where R is the conformally 
rescaled gravitational coupling, see the text) for the gravitino 
condensate field a in dynamically broken (conformal) Super¬ 
gravity models in the presence of a non-trivial de Sitter back¬ 
ground with cosmological constant A > 0 Q]. The Starobin- 
sky inflationary phase is associated with fluctuations of the 
condensate and gravitational field modes near the non-trivial 
minimum of the potential, where the condensate Uc 7 ^ 0 , and 
the potential assumes the value A > 0, consistent with super- 
symmetry breaking. The dashed green lines denote “forbid¬ 
den” areas of the condensate field values, violating the con¬ 
dition (fTTl). for which imaginary parts appear in the effective 
potential, thereby destabilizing the broken symmetry phase. 


loop cosmological constant A 0) : 


(19) 

with R denoting the fixed S'^ background we expand 
around {R = 4A, and the 4-dimensional Euclidean Vol¬ 
ume is 247r^/A^), and the a’s indicate the bosonic (gravi¬ 
ton) and fermionic (gravitino) quantum corrections at 
each order in A. 

The leading order term in A is then the effective action 
found in Q in the limit A —)> 0, 


Ea^o 

with 


247r^ 


Aq 


0 = 


A Z ( Aq 

Ai = - K - - 


On 


2471^ Ai 


( 20 ) 


( 21 ) 


^ The reader should notice that, upon the restriction US guaran¬ 
teeing the absence of imaginary parts in the one-loop effective ac¬ 
tion, the tree-level cosmological constant eg Ao < 0, while the 
one-loop one A > 0, as appropriate for a de Sitter background. 
Thus, Ao should not be confused with the current-epoch positive 
cosmological constant Aq, which we introduce later on, in section 
IV.21 when we discuss RVM (c/. II57I1 ). 


where 

Uq = (tI ^ 0.100 In ^ ^ +0-126^ , 


( 22 ) 


and 


Ctn 


_ ^4 /■ j;2 2Y 

— K (/ — cr^) 


1^0.027- 0.018 In 


(23) 


indicate the leading (as A —>• 0) contributions to the ef¬ 
fective potential from bosonic (graviton) and fermionic 
(gravitino) quantum fluctuations respectively, to one- 
loop order. Above, /r is a RG scale, associated with a 
short-distance proper time cutoff Q, not to be confused 
with the RG scale of the RVM fXc(t) (cf. Sect. II), which 
is such that the flow from Ultraviolet (UV) to Infrared 
(IR) corresponds to the direction of increasing /x, CTc de¬ 
notes the gravitino scalar condensate Oc (x k ip^) at 
the non-trivial minimum of the one-loop effective poten¬ 
tial {cf. Fig. [5]), K is the conformally-rescaled gravita¬ 
tional constant in the Jordan-frame SUGRA model of 
( 4 ^ . defined in (I14L corresponding to a non-trivial v.e.v. 
of the conformal (‘dilaton’) factor, (</)) 7 ^ 0 , assumed to 
be stabilized by means of an appropriate potential, lead¬ 
ing to the breaking of the conformal symmetry. In the 
case of standard A/" = 1 SUGRA, ((/>) = 0. 

The remaining (higher order in A) one-loop quan¬ 
tum corrections then, proportional to A and A^ may be 
identified respectively with Einstein-Hilbert i?-type and 
Starobinsky i?^-type terms in an effective action of the 
formS 


r ~ - 


2k2 


J' d^XyJ~g ^R — 2Ai^ -t- cx\ R -t- 0.2 R^ 


(24) 


where we have combined terms of order A^ into curvature 
scalar square terms. For general backgrounds such terms 
would correspond to invariants of the form Rf^upa , 

Rpv R^'' and R?, which for a de Sitter background all 
combine to yield R? terms 0. The coefficients ai and 


® The reader should recall at this stage that the sign of g in and 
the overall minus sign in front of the right-hand-side of is due 

to the Euclidean-signature formulation of the path integral and 
disappears upon analytic continuation back to the Minkowski 
space-time at the end of the computations, which is necessary in 
order to make contact with phenomenology/cosmology (see, e.g. 
©)• This should be understood in what follows, and especially 
in the context of linking the SUGRA model with the RVM in 
section lYl 

In the pure SUGRA case, with no dilaton frame functions, 
the fact that the Gauss Bonnet combination — 

AiRfiu R^^ -\- R^ is a total derivative in four space-time dimen¬ 
sion, implies that one can consider only the Ricci-scalar and 
Ricci-tensor squared terms as independent. This is not the case 
though in the conformal SUGRA case [S . 













di in (|24[) absorb the non-polynomial (logarithmic) in 
A contributions, so that we may then identify with 
CT via 

«i = Y(«f + «f), «2 = Y(ar+ar), (25) 


of dynamical breaking of supersymmetry In terms 
of the Starobinsky inflationary potential ([5]), the value 
A > 0 corresponds to the approximately constant value 
of this potential in the high (^-held regime (nip ^ 1) of 
Fig.m in the flat region where Starobinsky-type inflation 
takes place. Thus we may set 


where we note that ai is dimensionless whereas 02 has 
dimension of inverse mass squared. The coefficients 
z = 1, 2 in Eq. (|2^ can be computed using the results of 
derived via an asymptotic expansion: 


and 


of = 0.067 ~ 0.021 In + 


0.073 kV In 


= 0.029-b 0.014 In 


A 


(^) 


-0.0291n(-) , 


af = -0.083Ao -b 0.018 Aq In 


A 

V 


+ 


0.049 Ao In - 


3Ao 


af = 0.020-b 0.021 In ( 

\3/r2 


0.014 In - 


6 A 0 


(26) 


(27) 


A - 3 iJf (29) 

where Hj the (approximately) constant Hubble scale dur¬ 
ing inflation, which is constrained by the current data to 
satisfy @ or ([3]). In the SUGRA context under discus¬ 
sion Hj is linked to the scale of global SUSY breaking 
through Hi ~ f/Mpi. 

The effective Newton’s constant in (HI, after the im¬ 
position of (1351) . is then defined as 

= - , (30) 

Oil 

and from this, we can express the effective Starobinsky 
parameter ([5]) in terms of Keff as 

/3eff = — . (31) 

ai 

This condition thus makes a direct link between the ac¬ 
tion (1181) with a Starobinsky type action ([^ . Comparing 
with ([ 5 ]) we can determine the effective scalaron mass in 
this case: 


M = 


Stt ai 

3 02 


(32) 


To identify the conditions for phenomenologically accept¬ 
able Starobinsky inflation around the non-trivial minima 
of the broken SUGRA phase of our model, we impose 
first the cancellation of the “classical” Einstein-Hilbert 
space term R by the “cosmological constant” term Ai, 
i.e. that 

^ > 0 . (28) 

This condition should be understood as a necessary 
one characterizing our background in order to produce 
phenomenologically-acceptable Starobinsky inflation in 
the broken SUGRA phase following the first inflationary 
stage, as discussed in Q. This may naturally be under¬ 
stood as a TOneralization of the relation R = 2Ai = 0, 
imposed in Q as a self-consistency condition for the dy¬ 
namical generation of a gravitino mass in the flat (zero 
A) limit. 

From Eq. (l28ll it follows that the (positive) cosmolog¬ 
ical constant A > 0 satisfies the four-dimensional Ein¬ 
stein equations in the non-trivial minimum, and in fact 
coincides with the value of the one-loop effective poten¬ 
tial of the gravitino condensate at this minimum. As we 
discussed in this non-vanishing positive value of the 
effective potential is consistent with the generic features 


i? = 4A = 2Ai = -2k^ - 


An 


As we know, this mass parameter also sets the order of 
magnitude of the inflationary scale in the Starobinsky 
model. 

We may then determine the coefficients ai and 02 in 
order to evaluate the scale l/^/P ^ A4 of the effective 
Starobinsky potential given in Fig. [T] in this case, and 
thus the scale of the second inflationary phase. 

In 0 we searched numerically for points in the param¬ 
eter space such that: 


• The effective equations 




(33) 


are satisfied, together with the condition (l28ll . 

• the cosmological constant A is small and positive, 
satisfying (HU), and for phenomenological reasons 
it should be of order 


0 < A - 10"^°M^i , (34) 

to ensure the validity of our expansion in A, con¬ 
sistent with the phenomenology of Planck- satellite 
data 0. 






9 


• The scalaron mass should also be of order A4 ^ 
hence allowing us to achieve phe¬ 
nomenologically acceptable Starobinsky inflation in 
the massive gravitino phase, consistent with the 
Planck-satellite data 0. 

For K — K (he. for non-conformal supergravity), we 
were unable to find any solutions satisfying these con¬ 
straints. This of course may not be surprising, given 
the previously demonstrated non-phenomenological suit¬ 
ability of this simple model If we consider R, ^ k 
however, we find that we are able to satisfy the above 
constraints for a range of values El- 

In general, typical values obtained in phenomenological 
realistic conformal SUGRA models satisfy k ^ k {e.g. 
of order k/k = 0(10^ — 10^)), under the constraints (l29l) 
and in such a way that: 

A ~ ~ rny2 ~ <^c ~ j ^ /^. 

(35) 

Since the scale of SUSY breaking must be in the ball¬ 
park of the typical GUT scale associated to the inflation, 
namely \f] ^ 10 ^® GeV^ 10 “^ Mp;, from the above 
we have A ^ /Mpi ^ 10^^ GeV^. As a re¬ 

sult the scale of the gravitino is some two to three or¬ 
ders of magnitude below the GUT scale, that is to say, 
W 3/2 \/A > 10^^ GeV^ 10~^ Mpi. These values are 

compatible with both the combined Planck and Bicep2 
bound ([3]) and the typical mass of the gravitino in this 
framework Q . 

Exit from the inflationary phase is, of course, a com¬ 
plicated issue which we shall not discuss here at the level 
of the SUGRA model itself, aside from the observation 


® A comment concerning SUGRA models in the Jordan frame with 
such large values for their frame functions is in order here. In 
our approach, the dilaton 2<f> could be a genuine (dimension¬ 
less) dilation scalar field arising in the gravitational multiplet 
of string theory, whose low-energy limit may be identified with 
some form of SUGRA action. In our normalization the string 

( \-l/2 

coupling would begs=e'^ = lre/re) .In such a case, a value 

of it = K = 0(10®“^) K would imply a large negative v.e.v. 

of the (four-dimensional) dilaton field of order {<()) = —0(5) < 0, 
and thus a weak string coupling squared gs = 0(10“^), which 
may not be far from values attained in realistic phenomenolog¬ 
ical string models. On the other hand, in the Jordan-frame 
SUGRA models of [2, the frame function reads $ = = 

1-1(55 + HiHl) - lx ( - HS + Hi H- + h.c.), 
in the notation of 1461 for the various matter super fields of 
the next-to-minimal supersymmetric standard model that can 
be embedded in such supergravities. The quantity x is a con¬ 
stant parameter. At energy scales much lower than GUT, it is 
expected that the various fields take on subplanckian values, in 
which case the frame function is almost one, and hence re ~ re 
for such models today. To ensure re 3> re, and thus large val¬ 
ues of the frame function, 3> 1, as required in our analysis, 
one needs to invoke trasnplanckian values for some of the fields, 
and large values of x, which may indeed characterize the 
inflationary phase of such theories. A similar situation occurs for 
the values of the Higgs field (playing the role of the infiaton) in 
the non-supersymmetric Higgs inflation models [s^l . 


that it can be achieved by coherent oscillations of the 
gravitino condensate field around its minima and sub¬ 
sequent decays to radiation and matter fields (thus re¬ 
quiring detailed knowledge of the matter content of the 
SUGRA models in order to arrive at quantitative pre¬ 
dictions for the exit phase), or tunnelling processes a la 
Vilenkin [s^. However, in the next section we will show 
that the SUGRA model can be represented by an effec¬ 
tive running vacuum model along the lines indicated in 
Sect HU and from this point of view the exiting from the 
inflationary phase into the standard radiation phase can 
be guaranteed on very general grounds. 

Before doing so, though, we should make some impor¬ 
tant remarks concerning the presence of logarithms of the 
de Sitter scale A in the coefficients ai of the curvature 
terms of the effective action (l24ll . When one computes 
the effective action in a fixed de Sitter background, it 
is tempting to identify a A term with the Ricci scalar, 
which eventually will be allowed to depend on time. 
Thus, naively, the presence of logarithms would imply 
non-polynomial terms of the form i?lnii which would 
be problematic for any RVM interpretation of the exit 
from the inflationary phase, as it would contradict the 
spirit of the approach where only integer powers of the 
curvature terms would be allowed in the respective flow 
equations [l^, HH, EMI- Fortunately this is not the 
case. To understand this, we first remark that any effec¬ 
tive action obtained by integrating out massive degrees 
of freedom such as gravitino fields, that we restrict our¬ 
selves here, must consists for reasons of covariance and 
consistency of the weak gravitational fluctuations about 
the de Sitter background only of polynomial structures 
of the curvature tensors, for instance to fourth order in 
derivatives terms involving the squares of the Ricci scalar 
and Ricci tensors and covariant derivatives thereof. Any 
i?lni? term would be incompatible with the weak grav¬ 
ity perturbative expansion about a background, say of 
constant non-zero curvature. 

Thus, the coefficients ai and a 2 in the action ((Ml) are 
kept fixed, not undergoing temporal evolution, which is 
guaranteed by the fixing of the two free scales in the prob¬ 
lem n (1351) and A (1331) . Notice that the scale p, should 
not be confused with the subsequent RG scale that 
describes the cosmological evolution of the RVM vacuum 
(cf. Sect.ini)- Indeed, the scale /r first of all is a high 
energy cut-off. As already mentioned, it plays the role 
of a proper-time cutoff scale [^, appearing in the inte¬ 
gral representations of some C-functions that are part 
of the determinants arising in the path integral of the 
SUGRA action arising from integrating out massive spin 
3/2 (gravitino) and spin 2 (graviton) fluctuations about 
the de Sitter background. The scale fjL is therefore, in con¬ 
trast to Hc(t), an inverse renormalization group scale. Its 
value has to be fixed so as to guarantee SUGRA break¬ 
ing and to generate a fixed gravitino mass which should 
not depend on time. This implies that the spontaneous 
breaking of SUGRA and the inflationary phase are char¬ 
acterised by such fixed scales, which implies the time 
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independence of A (or, equivalently, the Hubble param¬ 
eter) during inflation, the gravitino mass, related to the 
gravitino condensate vacuum expectation value Uc, and 
thus the coefficients ai and 02 - On the other hand, inte¬ 
ger positive powers of A, appearing in the effective action 
may be replaced by higher order tensorial structures in¬ 
volving the square of the curvature tensors, which are 
allowed to vary with the cosmic time during the RVM 
phase after exit from inflation. Notice that microscopi¬ 
cally the exit phase is characterised by an unknown sort 
of phase transition, either through decays of the grav¬ 
itino condensates to matter parts and reheating of the 
inflated universe, or tunnelling, as mentioned previously, 
and thus using different RG running to relate various eras 
of the Universe after inflation is to be expected. 


V. “DECAY” OF EFFECTIVE VACUUM 
ENERGY: RUNNING VACUUM MODEL 

The main aim of this section is to demonstrate that 
there exists a family of time-dependent effective vacuum 
energy decaying models of running type, i.e. the class 
of the running vacuum models (RVM’s) introduced in 
Sect.im which characterize the evolution of the Universe 
from the exit of the Starobinsky inflationary phase till the 
present era. In fact, the RVM’s are able to interpolate on 
very general grounds the primeval de Sitter epoch with 
the late time de Sitter era, i.e. the dark energy one, where 
a much smaller cosmological constant essentially domi¬ 
nates. We shall follow the approach of the RVM outlined 
in Sect. m in which the vacuum energy density p\{H) 
varies with time through its dependence on H = H(t). 
The Hubble parameter, having dimension of energy in 
natural units, acts as the natural running scale via the 
RG equation Eq.Q. As mentioned in Sect. im only the 
even powers of H can be involved in that equation, ow¬ 
ing to the general covariance of the effective action. This 
is an important point to make possible a general QFT 
description of this RG approach and is essential for the 
connection with the SUGRA model under discussion. 

It is evident that the expansion quickly converges 
at low energies, where H is rather small - certainly much 
smaller than any particle mass. No other i?^"-term be¬ 
yond (not even can contribute significantly on 
the r.h.s. of equation (|1]) at any stage of the cosmologi¬ 
cal history below the GUT scale Mx ^ 10^® GeV, where 
presumably inflation occurs. 

On the other hand, if we want to deal with the physics 
of inflation and in general to the very early states of 
the cosmic evolution, we have to keep at least the term 
which in fact is the dominant term in the series o 
during the high energy regime. In contrast, the terms 
H^/Mf and above are less and less important because 
these higher and higher powers of H are suppressed by 
the inverse powers of the heavy fermion and boson masses 
in the GUT, as required by the Appelquist-Carazzone de¬ 
coupling theorem. Therefore, the dominant part of the 


series ® is expected to be naturally truncated at the 
term. Higher order terms should contain the bulk of the 
high energy contributions within Quantum Field Theory 
in curved spacetime, namely within a semi-classical de¬ 
scription of gravity near but (possibly a few orders) below 
the Planck scale. Models of inflation based on higher or¬ 
der terms inspired by the RG framework exist since long 
in the literature (see 0 ) as well as the unified inflation- 
dark energy framework of 0. For a more phenomeno¬ 
logical treatment unrelated to the RG, see [5J-[^. 


V.l. A distinct class of running vacuum models 


Based on the above arguments, it is natural to consider 
the case in which the highest power of the Hubble rate in 
the RG Eq. (|4l) is Integrating the RG equation pro¬ 
vides the simplest realization of RVM that can describe 
inflation and the various stages of the FLRW regime: 

Pa{H) = = A . (36) 


Here Cq is an integration constant (with dimension +2 
in natural units, i.e. energy squared) which can be fixed 
from the low energy data of the current universe [2l|, lU ■ 
On the other hand the dimensionless coefficients are given 
as follows: 


1 


E 


Mf 


(37) 


and 


a = 


1 H] 

9 ^ 


T. 


(38) 


At this point, we would like to make some comments 
which will hopefully make the reader appreciate the 
physical interpretation of the running vacuum scenario. 
The coefficient v behaves as the reduced (dimensionless) 
beta-function for the RG running of pA at low energies, 
whereas a plays a similar role at high energies. No¬ 
tice that the index i depends on whether bosons {B) 
or fermions {F) dominate in the loop contributions. Of 
course, since the coefficients {u, a) play the role of one- 
loop beta-functions (at the respective low and high en¬ 
ergy scales) they are expected to be naturally small be¬ 
cause Mf <C Mp[ for all the particles, even for the heavy 
fields of a typical GUT. Indeed, an estimate of v within a 
generic GUT is found in the range \v\ = 10“®— 10“^ 0. 
The dimensionless coefficient a is also small, |q:| 1, 

because the inflationary scale Hj is certainly below the 
Planck scale, see Eq. ([3]). From the observational view¬ 
point, utilizing a joint likelihood analysis of the recent 
supernovae type la data, the GMB shift parameter, and 
the Baryonic Acoustic Oscillations, it has been found 
\v\ = 0(10“®) [U, [ 23 , [13, 0, which is nicely in ac¬ 
cordance with the aforementioned theoretical expecta¬ 
tions as well as it insures a mild dynamical behavior of 
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the vacuum energy at low energies. As we have already 
stated in section im the Quantum-gravity corrections in 
the Starobinsky model have been found in the context of 
the RG analysis [H, through the appropriate beta- 
functions. The fact that the nature of the main coeffi¬ 
cients of both theories (running vacuum and Starobinsky) 
is based on the RG approach is a hint that perhaps there 
is a possible connection between the two models. 

Indeed as we will confirm below, this is the case. In 
particular, let us start with the effective action (I24L with 
coefficients (l25|) . (|26]) . (l27ll and the constraints (l28ll . (1^ . 
As we have seen, this action was obtained after integra¬ 
ting out both quantum-gravity (metric) fluctuations and 
massive gravitino fields. The action admits dynamical 
solutions of de-Sitter vacua with small cosmological con¬ 
stant A of order less than the GUT scale at early epochs 
of the Universe. At the non-trivial minimum, the effec¬ 
tive potential takes on a value of order A (c/. Fig. [2]). 
Around the minimum, we should replace A by an effec¬ 
tive Hubble parameter during inflation, Hj, Eq. (|2^ . 

At the end of the day the effective (dynamical) vacuum 
energy density, pa{H), during the inflationary phase of 
our SUGRA model can be extracted from the SUGRA 
effective action r (IMI), upon applying the constraint ( 1281 ) 
and analytically continuing the results back to Minkowski 
space-time signature. In particular, the effective poten¬ 
tial is defined as I4ff = —F —>■ / d'^Xy/^ pa{H). Doing 
so, we observe that the so obtained pA(iJ), remarkably, 
adopts precisely the generic RVM structure (1551) around 
that phase, in which the Ricci scalar - see Eq. dini) below 
— boils down to R ~ 12H^ since H remains (approxi¬ 
mately) constant in this phase. 

Some important remarks are in order here. The impo¬ 
sition of the constraint (1551) during the Starobinsky in¬ 
flationary phase implies, as already mentioned, that the 
correct phenomenology is attained as a result of the ef¬ 
fective gravitational coupling ([30]) that characterises that 
phase. So, if the constraint was an exact result, the effec¬ 
tive vacuum energy density of the SUGRA model would 
then correspond to the ^ 144 H'^ terms in (l24)) with 
(Uni) playing the role of the effective gravitational con¬ 
stant, 


SUGRA/^\exact _ *^2 „4 

Pa /constraint 2 ^ 

«eff 


18 af -b af 4 

r^F ^ > 

^eff -I- ai 


(39) 

where we used (1251) . (l26l) . (l27l) . The form (l39l) constitutes 
an admissible class of RVM (c/. Eq. (1551) 1. Notice that 
in Eq. (1391) there is no v term. This is important, in the 
sense that in such a model, as a result of the effective 
gravitational constant (1301) entering the game, which in 
this scenario [ 3 ] is viewed as the ^physicaV reduced Planck 
mass of order 10^® GeV, the gravitino mass and global 
SUSY breaking scales, (1551) . when expressed in terms of 
Keff are of order one, that is one encounters a Planck- 
scale gravitino. Despite this, the vanishing of v makes 
the renormalization-group equation (1391) a consistent one 
within the perturbative class of (1551) . 


However the above construction leads to the absence of 
a present-era (small, positive) cosmological constant cq. 
This arises from the fact that we imposed the constraint 
(1551) exactly. It may well be that such a condition leaves 
(non-perturbatively, when all the higher than one-loop 
contributions are taken into account) a very small (con¬ 
stant in cosmic time) contribution cq > 0 which is pre¬ 
served until the present day. Unfortunately, our one-loop 
construction does not allow us to explain the magnitude 
and the sign of this constant term, but this is equivalent 
to offering a solution to the cosmological constant prob¬ 
lem, which of course our approximate one-loop analysis 
cannot provide. While we do not have a quantitative 
calculation at this point, the above argument provides at 
least an interesting qualitative explanation, to wit: the 
origin of the current cosmological term p\ in the model 
might well be attributed to quantum (non-perturbative) 
effects in the SUGRA effective action, which prevent the 
complete cancellation (1551) from being realised. The con¬ 
stant residue cq is then transferred throughout the cosmic 
history and pops up in our days in the form of the tiny 
vacuum energy = (3/«:^)(co + uFIq) ~ 3co/k^. 

Under this assumption, then, the initial gravitational 
coupling K, and thus the Einstein term / d^x^—gR 
would enter the game during the exit phase from infla¬ 
tion 0. In such a case, in the exit phase, the effective 
vacuum energy of the SUGRA model at the inflationary 
phase should correspond to both ai and a 2 terms of (1241) . 
with the constraint (1281) failing by a tiny amount cq > 0 
corresponding to the present-era cosmological constant. 

^SUGRA = ^(^co + 6aiH^ + 72aa (40) 

where the explicit form of the coefficients ai, i = 1,2, 
given by Eqs. (EH), dMl) and (Ell), in which the scales p 
and A are fixed through (1551) and (1551) respectively. As 
mentioned already, fixing of the scale p and A, implies 
fixed values for the gravitino mass. 

It is worth stressing that the aforementioned ambigu¬ 
ity concerning the failure of the exact constraint (1551) can 
be avoided altogether by observing that the correspond¬ 
ing pa{H) ultimately derives from the RG equation ({4]) 
discussed in Sect. m It is therefore more appropriate, and 
elegant, if one performs the matching between the run¬ 
ning vacuum energies pa{H) in the SUGRA model and 
RVM by equating the corresponding “RG beta” func- 


The reader should bear in mind that, since during the infla¬ 
tionary phase the scalar degree of freedom of the Starobinsky 
action is slowly rolling, if there is inflation in the conformally 
rescaled metric m, there is also inflation in the initial metric. 
The Starobinsky inflation arguments are also not affected if a 
small contribution to the cosmological constant, of order of the 
present-era one, enters the effective action 0, as this is negligi¬ 
ble compared to the Hubble scale of inflation. 
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tions, 


d PA . 


———r = 6k ^ ai -ff ^ + 144 K a 2 H'^ ~ 
3 (af + af) + 18 (af + a^) H'^ 


(41) 


It is remarkable that the effective dynamical vacuum en¬ 
ergy density pa{H) associated to the SUGRA model un¬ 
der consideration turns out to follow the general RG of 
the RVM, see Eq. (|4]), in which the coefficients of the 
and H'^ terms can be computed precisely from the un¬ 
derlying SUGRA framework. 

From Eq. m it follows that pA = Pa{H) evolves 
(“runs”) with the (time-evolving) value of H. Such evo¬ 
lution will be studied in more detail in the next section, 
but is relatively small in the beginning, namely the vary¬ 
ing H is only slightly below the initial value given in 
Eq. (|29p . As a result the Universe can start with an ini¬ 
tial inflationary phase, which is dominated by the ^ 
term of m- However, well after the inflationary period 
the ~ iJ^-term takes over and remains in force until the 
present time, thereby providing a mild evolution of the 
current cosmological “constant”. 

Eventually one has to add to the effective action loop 
contributions from other matter fields, including parti¬ 
cle multiplicities, but at the moment we take a mass of 
order of the gravitino mass and shall comment on the 
possible additional effects below. The value of TO 3/2 is 
of order of the GUT scale, since it is proportional to 
the gravitino condensate through CTc ^ the 

latter being bounded from above by the GUT scale: 
y/oc < V7 ^ 10“^A/pi [see Eq. (IT71) ]. We must also keep 
in mind that for phenomenologically acceptable solutions 
of the broken SUGRA model the ratio r = k/k is forced 
to stay in the range r = 0(10^ — 10^). 

From the generic values (1351) we adopt here we find: 


dpA 

d\nH^ 


1.59 H'^ + 25.76 . 


(42) 


The integrated form of Ea. (H^ . yields of course the ef¬ 
fective vacuum energy density at the scale H, pa{H), 
Eq. (HUl) . within the current SUGRA scenario. We thus 
have: 


Pa{H) ~co + 1.59 -f 12.88, (43) 


where cq is the integration constant, which will play the 
role of the current-era cosmological constant, as already 
mentioned. The result naturally adopts the generic form 
of the canonical RVM, Eq. (l36l) with cq = k^co/ 3 and the 
effective values for the coefficients v and a given by 


l/eff ~ 0.53 ~ O 


^l/2\ 


Oeff ~ 4.30 ~ O 



(44) 

(45) 


Thus we observe that, within the context of the pure 
SUGRA model, where only the gravitino plays the role 


of “matter”, both coefficients are small, of typical or¬ 
der 10“®, in accordance with their interpretation as /3- 
function coefficients of the running vacuum energy den¬ 
sity. Let us also note that the above estimate for Ves 
nicely fits with the formal expression obtained in the 
different context of anomaly-induced inflation, where it 
also takes the structure dSZl), namely a quantity propor¬ 
tional to the (squared) ratio of a heavy particle scale 
(in gene ral a collection of them) to the Planck mass - 
see [lol for details. In the full SUGRA case, after its 
coupling to ordinary matter and radiation fields, other 
SUSY heavy fermions with masses near the SUSY break¬ 
ing scale should also contribute to r'es and a, and this 
should enhance the value of this coefficient very signifi¬ 
cantly, thus bringing the obtained result even closer to 
the situation studied in Ref. 0- Overall, such consider¬ 
ations in phenomenologically realistic SUGRA situations 
could bring these parameters to a range ^ 10“^ accessi¬ 
ble to current observations [U ■ 

We next remark that, in the general case where the 
parameters of the SUGRA model are varied from the 
generic values considered in (j35p . but within the allowed 
range, the values of and Oeff can also undergo some 
variation and the sign of t'es could change. However 
we stress that the sign of Oeff remains always positive, 
which is essential for a correct description of inflation. 
This can be seen explicitly by comparing the various log¬ 
arithms involved in the structure of the coefficient 02 of 
the LL^-term in Eq. (IdlTl . together with the size and sign 
of their respective numerical coefficients (c/. Eqs. (1251) . 
(Hi), (HZl)). The positivity of 02 is maintained through¬ 
out the physically allowed parameter space, and derives 
essentially from the fact that H] ^ a1 ^ ^ 3/2 ’ 

agreement with the generic result ([35]). 

Finally, let us note that the circumstance that VeS 
could have either sign can only affects the dynamics of 
the vacuum energy in the late universe. The phenomeno¬ 
logical implications for both signs have actually been ex¬ 
plored recently in [2l|, , see also 0 . 

The upshot of the above considerations points to the 
existence of a remarkable relation between the running 
vacuum model Ea. (l36l) with that of SUGRA Ea. (l43l) . In 
the next section we discuss the predicted inflationary sce¬ 
nario 00 in the context of the general RVM 0,0| , 
and provide some interesting phenomenology that can be 
tested for the low energy regime, namely for the current 
Universe. 


V.2. Running Vacuum Evolution: from current to 
inflationary era 

In this section we investigate the conditions under 
which the running vacuum model can provide an infla¬ 
tionary era. The point of this session is first to demon¬ 
strate that is, if one starts from an inflationary era, at 
an early epoch, obtained in the context of a microscopic 
model, such as the Starobinsky inflation induced in the 
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SUGRA model, then the RVM can smoothly connect it 
with the current era, characterised by a very small value 
of the vacuum energy, with a cosmology of ACDM type. 
We shall follow a bottom-up” approach, in which, by 
starting from a late epoch FLRW Universe and apply¬ 
ing RVM evolution (“backwards” in cosmic time, or, in a 
RG sense, an IR to UV flow), one arrives at an inflation¬ 
ary era in the early Universe. As we shall see, however, 
in this bottom-up approach, there is no unique way to 
identify the underlying microscopic model during the de 
Sitter era, which was to be expected in view of the rather 
generic features encapsulated in the RVM evolution. 

To this end, let us first reproduce the Friedmann equa¬ 
tions in the framework of a running pA. The resulting 
equations are expected to be formally equivalent to the 
AGDM case, inasmuch as the Gosmological Principle, 
which is embedded in the FLRW metric, perfectly al¬ 
lows the possibility of a time-evolving cosmological term. 
In general, the Einstein-Hilbert action is given by (here 
and in what follows we are back in Minkowski-signature 
space-time, described by a metric g^v)'- 


Sr,A = / d^Xy/^ 


2s;2 


{R - 2A) -f Cr, 


(46) 


where in our case pA{t) = A(t)/K^ represents the effec¬ 
tive vacuum energy, which is allowed to vary with the 
cosmic time (more specifically as a function of a dynam¬ 
ical cosmological variable that evolves with time), and 
£m is the Lagrangian of matter. Varying the action (l46ll 
with respect to the metric we arrive at 

Rui' - > (47) 


where the total is given by pR, with 

= -‘^dCm/dg^iv+g^iv £m the energy-momentum ten¬ 
sor corresponding to the matter Lagrangian. The extra 
piece is pA = A/k^, that is to say, the vacuum energy 
density associated to the presence of A(t) (with pressure 
PA = —Pa)- Let us remark that this equation of state 
(EoS) does not depend on whether the vacuum is dynam¬ 
ical or not. In contrast to other forms of dark energy, the 
vacuum is defined as that for which the EoS parameter 
is precisely w = — 1 in any circumstance. 

Modeling the expanding universe as a perfect fluid with 
velocity 4-vector field C/^, we obtain = p^ gfj,v + {Pm+ 
Pm) where pm is the density of matter-radiation 

and Pm = iXmPm is the corresponding pressure, in which 
uim is the EoS of matter. Obviously, takes the same 
form as with ptot = Pm + PA and ptot = Pm + PA = 
Pm PAt that is, = {jpm Pa) g^iv T (/^m 4” Pm)Rpd-dv- 

In the context of a spatially flat FLRW metric, we 
derive the Friedmann equations in the presence of a dy¬ 
namical A-term: 


K^Ptot = K^Pm + A = , (48) 


(49) 


and the Ricci scalar 

R = g^"Ru^ = + H), (50) 

where the overdot denotes derivative with respect to cos¬ 
mic time t. Note, that the Bianchi identities — Oj 

insure the covariance of the theory and, if the Newtonian 
coupling is strictly G =const., entail an energy exchange 
between vacuum and matter. 


Pm 4“ 3(1-1- UJrfi)R- Pm — PA • ( 51 ) 


Gombining equations (H5|) . (HHl) and m, we infer the 
basic differential equation that governs the dynamics of 
the Universe, namely the equation for the Hubble rate: 


H-\-^{l-\-uJm)H^ = ^ K^(l-fa;m)pA = (^ , (52) 

Inserting in it the expression (1361) for the dynamical vac¬ 
uum energy we arrive at the following equation: 


H 4- —(1 4- uJm)H‘^ 


, CO 


= 0 (53) 


The dynamics of this model has been thoroughly dis¬ 
cussed in [iB, Il3| , see also [l^ . We can summarize it as 
follows. First of all we identify the presence of an infla¬ 
tionary epoch (de Sitter phase) associated to the constant 
value solution = (1 — i/)Hf /a of Eq. (1551) . which is 
valid for the very early epoch of the universe (in which we 
can neglect cq/H^ <C 1). In this regime, solving Eg. (1551) 
we find 

H{a)=(^—^)^ (54) 

\ a J -pi’ 


where D is a positive constant of integration. For the 
early universe we assume that matter is essentially rela¬ 
tivistic, thus we take uJm = 1/3 at this point. Overall, 
one can see from (l54)) that for <C 1 the universe 

starts from an unstable inflationary phase [early de Sit¬ 
ter era, = (1 — v)Hj/a\ powered by the huge value 
iJ/ presumably connected to the scale of a Grand Uni¬ 
fied Theory (GUT). Well after the primeval inflationary 
era, specifically for ^ 1, the Universe enters 

the standard radiation phase. Subsequently the radia¬ 
tion component becomes subdominant and the matter 
dominated era appears. This is confirmed from the evo¬ 
lution of the vacuum energy and radiation energy densi¬ 
ties. If we neglect v and cojH^ in this early epoch, which 
is justified, we can insert (1551) into (1551) and we find: 


PA (a) = 


PI 


a [l + Da^Y 

Then solving (1511) we obtain: 

Pi Da^ 
Pr{a) = — 


a [l + DaY^ ' 


(55) 


K-'^Ptot = K^Pm - A = -2H - 37L^ 


( 56 ) 
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Here pi = ZH'j /is the critical density in the infla¬ 
tionary epoch. As it is obvious from the above expres¬ 
sions, there is no singularity in the initial state: the 
Universe starts at a = 0 with a huge vacuum energy 
density pi/a (and zero radiation) which is progressively 
converted into relativistic matter. In the asymptotic ra¬ 
diation regime we indeed retrieve the standard behavior 
Pr a~'^ with essentially negligible vacuum energy den¬ 
sity: PA ~ <C Pr- Graceful exit is, therefore, imple¬ 
mented. 

Subsequently the radiation component becomes sub¬ 
dominant and the matter dominated era appears. This 
is the point when the cq/H^ term in Ea. (l53l) surfaces and 
starts to dominate over aH^/Hf because the early de Sit¬ 
ter era is left well behind {H Hi). In this case Ea. (|36ll 
boils down to A(iJ) = Aq -I- — Hq) which corrects 

the concordance ACDM model a posteriori. Notice, that 

Ao = 3co -f 3uH^ (57) 


at low energies this implies that only the dynamical part 
^ is active and may lead to interesting phenomeno¬ 
logical irnplications for the dynamical DE of the current 
universe [2l|, [H, 1^1 


V.3. Geometrical description: RVM versus 
Starobinsky 

Finally, let us focus now on some aspects of the infla¬ 
tionary era that are especially relevant for the present 
study. As we have already mentioned, in this epoch we 
have a de Sitter solution ~ (1—z/)iJ|/a =const. Now, 
as previously indicated, ~ 0 in this period and hence 
R ~ 12iJ^. Finally, neglecting the matter component 
from the action (gSl), which is justified in the inflation¬ 
ary period, and using A(i7) ~ 3aH'^/H‘f [see Eg. (1551) ] we 
schematically find 


is the vacuum (cosmological constant) energy density at 
the present time, which is positive, and should not be 
confused with the negative tree-level Aq of the SUGRA 
model ®. This can be understood by studying the 
evolution of the universe at a time after recombination, 
therefore consisting of dust (w^ = 0) plus the running 
vacuum fluid with H <C Hj. In this case, using d/dt = 
aHd/da, we can rewrite Ea. (l53l) as 

rfM^ 

a — -1- 3(1 — iy)H^ — 3 cq = 0 , 

da 

The solution satisfying the boundary condition H — Hq 
at present (a = 1) is: 


H^{a) = 


Hi 

\-vV 




Note that the aforementioned boundary condition fixes 
the value of the parameter cq as follows: cq = Hq{VI/,^ — v). 
For V — 0 we correctly recover the behavior of the 
AGDM. However, for small v the Universe possesses a 
mildly evolving vacuum energy that could appear as dy¬ 
namical dark energy without invoking spurious scalar 
fields. Furthermore, the above vacuum model is in agree¬ 
ment with the latest cosmological data and it predicts a 
growth rate of clustering which is in agreement with the 
observations (for more details concerning the late dynam¬ 
ics, see [U [2J, ang). 

Let us note that the main stage of the cosmic evolu¬ 
tion where we can match the SUGRA model of Sect. IV 
with the RVM is the early period comprising inflation 
and the incipient radiation epoch, to which it leads af¬ 
ter graceful exit, as described in this section. Later on 
the microscopic description of the SUGRA model is more 
difficult to analyze and we adopt here the point of view 
that the subsequent effective behavior of the Universe 
still follows the RVM flow dictated by the general RG 
equation (|3]), which, to order H'^, entails the dynamical 
vacuum energy density (1551) . As previously mentioned. 


Sr,k = J d'^x^/^ 


^R-Pa{H) 


H^ 


— I dxV^g^R-6aj^ 


(58) 


This demonstrates our point that an inflationary vacuum 
can be connected smoothly, under the RVM, with a late 
epoch ACDM Universe. However, there is no unique way 
by means of which we can associate the inflationary era 
RVM effective action (1551) to a microscopic model, which, 
as already mentioned, is to be expected due to the generic 
features of the RVM that describe classes of models and 
therefore may correspond to more than one microscopic 
theories, as far as the exit from inflationary phase is con¬ 
cerned. 

An interesting point concerns Eq. (1551) if one replaces 
H'^ by the square of the Ricci scalar. In this case one may 
write Sr^a ~ f d'^x^/^ ^ (^R — ’ Notice that, 

since a > 0 in our case, the RVM model is not formally 
and directly equivalent to a Starobinski-type model, for 
which the effective Lagrangian has the form ([5]) corre¬ 
sponding to a negative a coefficient in ((55)) . This point 
has also been discussed in [I^ . The root of the problem 
lies in the fact that the metric tensors of the two mod¬ 
els, and (1551) are different, related by a non-trivial 
conformal transformation m involving the linearising 
Hubbard-Stratonovich field p, which plays the role of the 
“physical” inflaton. The RVM metric is identified with 
the Einstein-frame metric in while the original 
one-loop effective SUGRA action is described in terms of 
the Pfii, metric. Nevertheless, contact with Starobinsky- 
type models, like the one induced within the context of 
SUGRA model examined here, can be achieved by ob¬ 
serving that it is precisely the passage from the Einstein 
to Jordan-frame actions, via ®, which guarantees the 
opposite sign, relative to the Ricci scalar term, of the 
effective potential ([5]) of the Hubbard-Strstonovich infla¬ 
ton field in ®. Upon making the identification for large 
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2> 1 

rr4 (l- 

® 

where in the SUGRA model the scalaron mass scale 
M is given by (|5^ . one obtains the connection of the 
RVM with the microscopic Starobinsky-type inflationary 
SUGRA model. The exit from the inflationary phase, 
then, which in fig. [T] corresponds to the region of small 
K (/? < 4, which in the context of the SUGRA model would 
require detailed knowledge of the matter content of the 
theory, is then “effectively” described by the RVM evo¬ 
lution with the initial condition (I59L that “fills up” the 
missing details in the exit-phase of the evolution in a 
rather generic manner. Below we compare the two cos¬ 
mological models at the dynamical level. 



FIG. 3. The RVM effective potential an^U/H] (solid line) 
versns the scalaron field nup. In order to produce the curve 
we utilize A4 ~ Mx ~ 10^®Gev, Hj ~ 0.81 x 10^‘^Gev [see 
Eq.Q] and a ~ 10~^. The dashed line corresponds to the 
Starobinsky effective potential (see Fig. [!}. 


V.4. Scalar field description: RVM versus 
Starobinsky 

Although the fundamental origin of the RVM has a 
root in the general structure of the effective action of 
QFT in curved space-time, we cannot provide the latter 
at this point, see [l^ for an explanation. However, we 
can resemble it via an effective scalar field <j) using a field 
theoretical language. We may call this scalar field (j) as 
vacuumon. Based on Friedmann’s Eas. (l48ll - (l4^ and fol¬ 
lowing standard lines, namely ptot = P 4 > = ^ + G((^) 

and ptot = P 4 , = /‘2‘ — U (0) we arrive at 


(?!)(o) = — sinh ^ Y/Da^ , 

= iln (VDo^ + \/Da^ -h l) . (64) 

In this context, utilizing Eas. (l6T]l - (l63l) the effective po¬ 
tential is given by 


[/(«) = 


iJ? 3-1- Da^ 






2 • 


which implies 



(60) 



tt2 

uw = 

1 

^ aH'\ 

(61) or 

1 

/+3H ■ 


{1 + Da‘^Y' 


3 + sinh^(K(/)) 


where U{(j)) is the effective potential, H = aHH and 
prime here denotes derivative with respect to the scale 
factor. Integrating Eg. (16(1 we have 


(62) 


Now, for Wm = 1/3 the Hubble parameter (iMl) takes 
the form 




C0sh^(K(/)) 


cosh'*(K(/)) 


(65) 


( 66 ) 


(67) 



At this point we would like to pose the following dynam¬ 
ical question: Under what conditions the RVM potential 
is equal to that of Starobinsky, namely U{4>) = t4ff(p) in 
Eq. 

Equating the right-hand-side of Eqs.([5]) and (1571) . after 
some calculations, we can express the vacuumon field in 
terms of the scalaron 




iy/2 

a) ^Da^ + 1 


Hi 


(63) 


4>{p) = -In X(V’) + Vxip)'^ - 1 
K L 


( 68 ) 


Notice that, we have set ^ = 0 in Ea. (l54l) . which is not 
important for the study of the early universe. Inserting 
Eg. (1631) into Eg. (1621) and performing the integration in 
the interval [0, a] we find 


where xif) > 1 as easily shown (hence no restrictions on 
the scalar fields) and it is given by 


xip) = 


1 -b V'l -b 8F{p) 
2F{p) 


1/2 


(69) 
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F{(p) = 


3aM^ (^1-e-Vl^‘^y 

4fff 


aK 

Hf 


Vesiv) > 0 . 


In Fig. [3] we present the RVM effective potential 
aK^U/Hj (solid curve) as a function of K(p. In the same 
figure we plot the effective Starobinsky (dashed curve) 
potential which is showed in Fig.[TJ From the comparison 
it becomes clear that, although the RVM and Starobinsky 
models live in different geometrical backgrounds, namely 
GR and R^, the two models are similar from the point of 
view of those features of inflation that can be described 
by an effective scalar-field dynamics. However, in other 
important aspects they are different. We should mention 
that the RVM model provides a simple description of the 
graceful exit and reheating problem, see ©[BE Hi 
for details. As for the Starobinsky model, the reheating 
of the universe after the exit of the inflationary phase it 
has been discussed for example in [B| • In contrast to the 
Starobinsky model, a general effective action from where 
the RVM can be derived is not known [B| , and currently 
this has been achieved only in some cases M- 


VI. CONCLUSIONS 

In the light of the latest Planck-|-Bicep2 results @ it 
has been proposed that the Starobinsky inflation plays a 
key role because it fits quite well the Cosmic Microwave 
Background (CMB) data on inflation. In the present pa¬ 
per we have further investigated the class of the running 
vacuum models (RVM) |14j (based on renormalization- 
group approach in curved spacetimes) and their impli¬ 
cations on the inflationary universe [lO . m. In partic¬ 
ular, we have addressed the possibility that they can 
mimic both the original Starobinsky model and the spon¬ 
taneously broken SUGRA models based on dynamically 
induced gravitino condensates Q . 

We have shown that the vacuum energy density pk{H) 
of these SUGRA models can be expressed as an even 
power series ([4]) of the Hubble parameter, which can be 
naturally truncated at the term. This is exactly the 
generic form expected in the simplest class of running 
vacuum models and therefore we can ap ply the known im¬ 
plications of these models for inflation |l6l. [II- Namely, 
after computing the modified form of the Friedmann 
equation, we find that the physics of inflation (which 
in our case occurs for H ~ Hj, a value associated to 
the spontaneously broken SUGRA model) is mainly de¬ 
scribed by the H'^-term. Furthermore, being of order 
R^, we can trace some relationship of this model with 
Starobinsky inflation, although of course there is not a 
full identification or equivalence. Most noticeably we 
point out the distinguishing feature that within the en¬ 
tire class of running vacuum models - and hence, in par¬ 
ticular, the SUGRA model that we have studied (which 
adapts to the same pattern) - the RVM performs suc¬ 
cessful graceful exit from the inflationary phase into the 


standard radiation regime El IB- This feature is charac¬ 
teristic of the running type of vacuum models, in contrast 
to the original Starobinsky model. Nevertheless, we have 
also shown here that the RVM model admits a scalar 
field description as well, via the vacuumon field, and its 
potential can be made equivalent to the Starobinsky po¬ 
tential upon appropriate scalar field redefinitions, despite 
the fact that the geometric backgrounds of the two mod¬ 
els are very different. This dynamical equivalence implies 
that the two models should provide the same inflationary 
features, at least in all the aspects that can be described 
through an effective scalar field potential. Not so in other 
aspects which may differ from one model to the other. In 
particular let us emphasize that the Starobinsky model 
derives from a local effective action whereas the struc¬ 
ture of the effective action in the general case RVM is 
not presently known, except in particular cases in which 
it is found to be non-local. 

The low energy physics, on the other hand, and in 
particular the evolution of the Universe in the current 
epoch, is determined by the constant additive term of 
Pa{H) and the power , which provides a remnant dy¬ 
namical evolution still in our days, which is of the form 
Pa{H) — Pa + (3 — Hq). Such evolution is mild 

because the coefficient of is small (it is interpreted 
as the /3-function coefficient of the running vacuum en¬ 
ergy at present). The signature of the RVM at present 
is precisely that mild quadratic dynamical behavior of 
the vacuum energy density around the current value p^ 
which is parameterized by the small parameter v. The 
model has been thoroughly put to the test recently and 
it allows values of \v\ — 0(10“^) [Bllljlilillll)- On the 
other hand, its successful performance in describing the 
physics of the early Universe (in particular the graceful 
exit of the inflationary phase into the standard radiation 
one) is also quite encouraging, especially after realizing 
that specific QFT models lead to this kind of behavior. 
In this paper we have shown that SUGRA models with a 
dynamically induced massive gravitino phase lead to the 
RVM behavior and therefore provide a strong support for 
a fundamental description of the cosmic history. 

Finally, we would like to stress that, in the context of 
the running vacuum model, the universe evolution, and 
especially its accelerated phase either during inflation or 
at late times, is not attributed to an ad hoc scalar field, 
or to a modification of the gravitational interaction, but 
rather arises from the modification of the vacuum itself, 
which is endowed with a dynamical nature. Remarkably, 
the SUGRA framework studied here provides a concrete 
realization of this possibility within the fundamental con¬ 
text of quantum field theory in curved spacetime. 

Acknowledgments. SB acknowledges support by 
the Research Center for Astronomy of the Academy 
of Athens in the context of the program ‘‘‘‘Tracing the 
Cosmic Acceleration". The work of N.E.M. is sup¬ 
ported in part by the London Centre for Terauniverse 
Studies (LOTS), using funding from the European Re- 



17 


search Council via the Advanced Investigator Grant 
267352 and by STFC (UK) under the research grant 
ST/L000326/1. The work of JS has been partially sup¬ 


ported by FPA2013-46570 (MICINN), Consolider grant 
CSD2007-00042 (CPAN), 2014-SGR-104 (Generalitat de 
Catalunya) and MDM-2014-0369 (ICCUB). 


[1] For an up-to-date review of inflationary models, and com¬ 
parison to recent data, see: J. Martin, C. Ringeval and 
V. Vennin, Phys.Dark Univ. 5-6 (2014) 75. 

[2] For Planck constraints on inflationary models, 
see P. A. R. Ade et al. [Planck Collaboration], 
arXiv:1303.5082 [astro-ph.COj; for a general survey 
of Planck results, including measurement of cosmological 
parameters, see P. A. R. Ade et al. [Planck Collabora¬ 
tion], arXiv: 1303.5062 [astro-ph.CO]. P. A. R. Ade et al. 
[Planck Collaboration], Astron. Astrophys. 571 (2014) 
A16. 

[3] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980). 

[4] P. A. R. Ade et al. [BICEP2 Collaboration], 
arXiv:1403.3985 [astro-ph.CO]. 

[5] BICEP2 and Planck Collaborations (P.A.R. Ade (Cardiff 
U.) et al.). Joint Analysis of BICEP2/Keck Array and 
Planck Data Phys.Rev.Lett. 114 (2015) 10, 101301. 

[6] J. Alexandre, N. Houston and N. E. Mavromatos, Phys. 
Rev. D 88, 125017 (2013). Int. J. Mod. Phys. D 24 (2015) 
1541004. 

[7] J. Alexandre, N. Houston and N. E. Mavromatos, Phys. 
Rev. D 89 (2014) 027703. 

[8] J. Ellis and N. E. Mavromatos, Phys. Rev. D 88 (2013) 
085029. 

[9] 1. L. Shapiro and J. Sola, Phys. Lett. B530 (2002) 10; 
Russ. Phys. J. 45 (2002) 727, Izv. Vuz. Fiz. 2002N7 
(2002) 75. 

[10] J. Sola, J. Phys. Math. Theor. A41 (2008) 164066. 

[11] 1. L. Shapiro and J. Sola, JHEP 0202 (2002) 006; 
Phys.Lett. B475 (2000) 236. 

[12] 1. L. Shapiro and J. Sola, Phys. Lett. B 682, 105 (2009); 
e-Print: arXiv:0808.0315 

[13] J. Sola, H. Stefancic, Phys. Lett. B624 (2005) 147; Mod. 
Phys. Lett. A21 (2006) 479; 1. L. Shapiro, J. Sola, H. 
Stefancic JCAP 0501 (2005) 012. 

[14] J. Sola, J. Phys. Conf. Ser. 453 (2013) 012015 
arXiv: 1306.1527 ; AIP Conf.Proc. 1606 (2014) 19 
arXiv: 1402.7049). 

[15] J. Sola, and A. Gomez-Valent, Int. J. of Mod. Phys. D24 
(2015) 1541003. 

[16] J. A. S. Lima, S. Basilakos and J. Sola, Mon. Not. Roy. 
Astron. Soc. 431, 923 (2013). 

[17] S. Basilakos, J. A. S. Lima and J. Sola, Int. J. Mod. Phys. 
D 22, 1342008 (2013); Int.J.Mod.Phys. D23 (2014) 
1442011. 

[18] E. L. D. Perico, J. A. S. Lima, S. Basilakos and J. Sola, 
Phys. Rev. D 88, 063531 (2013). 

[19] J. Sola, Int. J. Mod. Phys. D24 (2015) 1544027. 

[20] J.A.S. Lima, S, Basilakos, J. Sola, Gen. Rel. Grav. 47 
(2015) 40. 

[21] A. Gomez-Valent, J. Sola, S. Basilakos, JCAP 01 (2015) 
004. 

[22] A. Gomez-Valent, J. Sola, Mon.Not.Roy.Astron.Soc. 448 
(2015) 2810. 

[23] A. Gomez-Valent, E. Karimkhani, J. Sola, JCAP 12 
(2015) 048. 


[24] S. Basilakos, M. Plionis, and J. Sola, Phys. Rev. D80 
(2009) 083511. 

[25] J. Grande, J. Sola, S. Basilakos and M. Plionis, JGAP 
08, 007 (2011). 

[26] H. Fritzsch, J. Sola, Class. Quant. Grav. 29 (2012) 
215002; J. Sola, Int. J. Mod. Phys. A29 (2014) 1444016. 

[27] S. Basilakos, D. Polarski, and J. Sola, Phys. Rev. D86 
(043010) 2012. 

[28] S. Basilakos and J. Sola, Mon. Not. Roy. Astron. Soc. 
437 (2014) 3331; Phys.Rev. D90 (2014) 023008. 

[29] G. Espana-Bonet et al., JCAP 0402 (2004) 006; 
Phys.Lett. B574 (2003) 149. 

[30] J. Sola, Running Vacuum in the Universe: current phe¬ 
nomenological status, to appear in the Proc. of the 
14th Marcel Grossman Meeting, Rome (2015), e-Print: 
arXiv:1601.01668 

[31] J. Sola, A. Gomez-Valent, J. de Cruz Perez, Astrophys. 
J. 811 (2015) L14. 

[32] J. Sola, et ah, e-Print:arXiv: 1602.02103, 1606.00450. 

[33] J. R. Ellis, N. E. Mavromatos and D. V. Nanopoulos, 
Gen. Rel. Grav. 32, 943 (2000) gr-qc/9810086 ; E. Gra- 
vanis and N. E. Mavromatos, Phys. Lett. B 547, 117 
(2002) hep-th/0205298 ; New J. Phys. 6, 171 (2004) 
gr-qc/0407089 . 

[34] G.W. Misner, K.S. Thorn, J.A. Wheeler, Gravitation 
(Freeman, San Francisco, 1973). 

[35] P. G. W. Davies, S. A. Fulling, S. M. Christensen and 
T. S. Bunch, Annals Phys. 109, 108 (1977); T. S. Bunch 
and P. C. W. Davies, Proc. Roy. Soc. Lond. A 360, 
117 (1978); N. Birrell and P. Davies, Quantum Fields in 
Curved Space, (Cambridge Monogr.Math.Phys., 1982). 

[36] A. Vilenkin, Phys. Rev. D 32, 2511 (1985). 

[37] K. S. Stelle, Gen. Rel. Grav. 9, 353 (1978); B. Whitt, 
Phys. Lett. B 145, 176 (1984). 

[38] E. J. Copeland, C. Rahmede and 1. D. Saltas, Phys.Rev. 
D91 (2015) 103530. 

[39] For reviews see: D. F. Litim, PoS QG-Ph (2007) 024; 
M. Reuter and F. Saueressig, New J. Phys. 14, 055022 
( 2012 ). 

[40] S. Weinberg, (1979), in General Relativity: An Einstein 
centenary survey (ed. S. W. Hawking and W. Israel, 
1979), 790- 831. 

[41] J. Ellis, D. V. Nanopoulos and K. A. Olive, Phys. Rev. 
Lett. Ill (2013) 111301; JCAP 1310 (2013) 009. 

[42] K. Nakayama, F. Takahashi and T. T. Yanagida, 
Phys.Lett. B725 (2013) 111; JCAP 1308 (2013) 038; 
R. Kallosh and A. Linde, JCAP 1306, 028 (2013); 
W. Buchmuller, V. Domcke and K. Kamada, Phys.Lett. 
B726 (2013) 467; M. A. G. Garcia and K. A. Olive, JCAP 
1309 (2013) 007; F. Farakos, A. Kehagias and A. Riotto, 
Nucl. Phys. B 876, 187 (2013); D. Roest, M. Scalisi and 
1. Zavala, JCAP 1311 (2013) 007; S. Ferrara, R. Kallosh, 
A. Linde and M. Porrati, Phys.Rev. D88 (2013) 085038. 

[43] P. Van Nieuwenhuizen, Phys. Rept. 68 (1981) 189; 
D. Z. Freedman and A. Van Proeyen, Cambridge, UK: 
Cambridge Univ. Pr. (2012) 607 p; 




18 


[44] J. Jaeckel and C. Wetterich, Phys. Rev. D 68 (2003) 
025020 hep-ph/0207094 . 

[45] S. Ferrara, R. Kallosh, A. Linde, A. Marrani and A. Van 
Proeyen, Phys. Rev. D 82, 045003 (2010); ibid. D 83, 
025008 (2011). 

[46] M. B. Einhorn and D. R. T. Jones, JHEP 1003 (2010) 
026. 

[47] S. Deser and B. Zumino, Phys. Rev. Lett. 38 (1977) 1433. 

[48] Z. Komargodski and N. Seiberg, JHEP 0909 (2009) 066. 

[49] 1. L. Buchbinder and S. D. Odintsov, Class. Quant. Grav. 

6 (1989) 1955; see also S. D. Odintsov, Phys. Lett. B 213, 

7 (1988). 

]50] E. S. Fradkin and A. A. Tseytlin, Nucl. Phys. B 234 
(1984) 509. 

]51] E. S. Fradkin and A. A. Tseytlin, Nucl. Phys. B 234 
(1984) 472. 

[52] E. Witten, Nucl. Phys. B 188, 513 (1981). 

[53] F. L. Bezrukov and M. Shaposhnikov, Phys. Lett. B 659 
(2008) 703; F. Bezrukov, A. Magnin, M. Shaposhnikov 


and S. Sibiryakov, JHEP 1101 (2011) 016. 

[54] J. A. S. Lima and J. M. F. Maia, Phys. Rev. D 49, 5597 
(1994). 

[55] J. A. S. Lima and M. Trodden, Phys. Rev. D 53, 4280 
(1996), astro-ph/9508049 

[56] J. M. F. Maia and J. A. S. Lima (2000), unpublished. See 
also, J. M. F. Maia, Some Applications of Scalar Fields in 
Cosmoloqy, PhD thesis (in Portuguese), Sao Paulo Uni¬ 
versity (2000). 

[57] S. Carneiro and R. Tavakol, Gen. Rel. Grav. 41, 2287 
(2009). 

[58] J.A.S. Lima, S. Basilakos, J. Sola, Eur. Phys. J. CTO 
(2016). 

[59] D.S. Gorbunov, A.G. Panin, Phys.Lett. B700 (2011) 157; 
E.V. Arbuzova, A.D. Dolgov, L. Reverberi, JCAP 1202 
(2012) 049; H. Motohashi, A. Nishizawa, Phys.Rev. D86 
(2012) 083514. 


